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ABSTRACT 

We present the results of a deep search for pulsations from the nearby isolated neutron star RX J 1856.5-3754 
using the 450ks Director's Discretionary Time Chandra observation completed on 2001 Oct 15. No pulsations 
were detected. We find a 99% confidence upper limit on the pulsed fraction of ~ 4.5% for worst-case sinusoidal 
pulsations with frequency < 50 Hz and frequency derivatives -5 x 10~ 10 < / < Hz s" 1 . The non-detection of 
pulsations is most likely due to an unfavorable viewing angle or emitting geometry. Such emitting geometries are 
much more likely to occur for more "compact" neutron stars which show increased gravitational light-bending 



effects. In this case, the non-detection implies a radius/mass ratio for RX J 1856.5-3754 of R/M < lOkmM 
Subject headings: stars: individual (RX J 1856.5-3754) — stars: neutron — X-rays: stars 
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tions and placed an upper limit on the pulsed fraction (for fre- 
quencies < 40 Hz) of ~ 8%. Similar searches have been un- 
dertaken with archival ROSAT and ASCA data by Pons et al. 
(2002) who placed a 50% confidence limit on the pulsed frac- 
tion of ~ 6%. It should be noted that three other radio quiet 
isolated NS candidates have been reported to show X-ray pul- 
sations: RX J0720.4-3125 with a period of 8.4 s and a pulsed 
fraction of - 12% (Haberl et al. 1997), RX J0420.0-5022 with 
a period of 22.7 s and a pulsed fraction of ~ 40% (Haberl et al. 
1999), and RBS 1223 with a period of 5.16 s and a pulsed frac- 
tion of ~ 20% (Hambaryan et al. 2002). 

In 2001 October, RXJ1856 was observed with the Chan- 
dra LETGS for an additional 450 ks of Director's Discretionary 
Time (DDT) over the course of ~ 7.25 days. The purpose of the 
observation was to search for features in the spectrum — in par- 
ticular, due to heavy elements in the NS atmosphere — in order 
to constrain the NS equation of state (EOS). However, these 
data are of sufficient time resolution to allow a deep search for 
pulsations as well. The data were made public soon after the 
observation was completed. We conducted brute-force period 
folding and advanced Fourier analyses of this data in order to 
try to identify the spin parameters of the neutron star. In addi- 
tion, we have analyzed archival observations of RXJ1856 from 
both ROSAT and Chandra in order to specify a consistent set of 
upper limits to the pulsed fraction. 

2. OBSERVATIONS AND DATA PREPARATION 

The DDT observation comprised Obslds 3382, 3380, 3381, 
and 3399 in chronological order beginning on 2001 October 8 
and ending on 2001 October 15 (see Table 1) and were taken 
in the standard LETGS configuration. The total exposure time 
was 450 ks with significant gaps between the various Obslds 
giving a total duration of ~ 626 ks. 

The data were prepared in a very simple manner. We ex- 
tracted all events from within ~ 2" of the centroid of the zero- 
order image from each Obsld. The arrival times of these events 
were transformed to the Solar System Barycenter using the 
standard Chandra Interactive Analysis of Observations ( CIAO, 
v2.2) tool axBary, the nominal position of RXJ1856 from 
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1. INTRODUCTION 

The ROSAT observatory provided us with at least seven soft 
X-ray emitting point sources which are likely to be radio-quiet 
thermally-emitting isolated neutron stars (NSs; for reviews see 
Treves et al. 2000; Motch 2002). The nearest of these systems 
is RX J1856.5-3754 (hereafter RXJ1856), which was first ob- 
served as an unidentified point source in the EINSTEIN Slew 
Survey (Elvis et al. 1992) and later identified as an isolated NS 
by Walter, Wolk, & Neuhauser (1996). Subsequent observa- 
tions in the optical and ultraviolet (UV) have strengthened this 
claim by detecting a faint blue optical counterpart (Walter & 
Matthews 1997), measuring its proper motion (Walter 2001), 
and identifying a bow-shock H a nebula around the source (van 
Kerkwijk & Kulkarni 2001). Walter (2001) also used Hubble 
Space Telescope (HST) data to measure a parallax distance of 
61^pc to RXJ1856. However, Kaplan, van Kerkwijk, & An- 
derson (2002) have recently re-analyzed the same data and de- 
termined a distance of 140 ±40pc — the reasons for the large 
discrepancy are, as of yet, unclear. 

Despite these measurements, the nature of RXJ1856 remains 
a mystery. Is it young or old? Is it an ordinary radio pulsar 
beaming away from us, a low magnetic field neutron star ac- 
creting from the interstellar medium, or a high-field "magne- 
tar"? A key to answering these questions lies in detecting pul- 
sations from this source, from which we could then determine 
its spin-period, age, braking torque, and magnetic field strength. 
If an isolated neutron star is sufficiently magnetized, the result- 
ing anisotropies in the temperature distribution of the crust can 
produce hotspots of X-ray emission on the surface of the star. 
If the magnetic axis is mis-aligned with the rotation axis, these 
hotspots may produce a periodic modulation of the soft X-ray 
intensity at the rotation period of the NS. Even a non-detection 
of these pulsations can provide useful constraints on the mass, 
radius, and emitting geometry of the NS. 

Searches for these pulsations were carried out by Bur- 
witz et al. (2001) using a ~ 56ks Chandra X-Ray Observa- 
tory observation of RXJ1856 taken on 2000 March 10 using 
the Low Energy Transmission Grating Spectograph (LETGS 4 , 
Brinkman et al. 2000). Burwitz et al. (2001) found no pulsa- 
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Walter (2001) of RA(J2000) = 18 h 56 m 35 ?5 and DEC(J2000) = 
-37° 54' 36"8, and the preliminary Level 2 orbital ephemeris 
file distributed with the data release. We binned the 90134 re- 
sulting events into 1 .5 ms time bins to create a 420 million point 
time series that was used in the coherent Fourier analysis de- 
scribed in §3.1. 

A wiring error in the HRC causes each event to be tagged 
with the arrival time of the previous event. Unfortunately not 
every event is telemetered to the ground, so the arrival times of 
recorded events are typically in error by a few ms 5 . In order 
to improve the accuracy of these arrival times and to restrict 
the number of counts so as to make a brute-force period fold- 
ing search computationally feasible, we also reprocessed the 
data using the time filtering scheme described by Tennant et al. 
(2001). By shifting all arrival times in the Level 1 event files 
back to the previous event, filtering the data spatially, barycen- 
tering, and then keeping only those events that arrived within 
1 ms of the previous event, we are guaranteed to have all ar- 
rival times accurate to < 1 ms at the cost of a loss of ~ 95% of 
the events that were used in the Fourier analysis (see Table 1). 
These time filtered events were used in the folding search de- 
scribed in §3.2. 

In order to verify preparation methods for the data, we pro- 
cessed events from the LETGS observation of the Crab pulsar 
(Obsld 759) in the same manner as described above. The Crab 
pulsar was easily detected in the data at the expected barycen- 
tric rotation period and with the expected pulse profile given the 
HRC wiring error. In addition, significant detection of all spin 
harmonics up to at least the eighth (<~ 239 Hz), indicate that 
fast periodicities (> 100 Hz) are detectable — albeit drastically 
suppressed — in unfiltered HRC data given enough signal-to- 
noise. 

3. DATA ANALYSIS 

In order to maximize our sensitivity to coherent pulsations 
with a variety of pulse shapes, we searched the data using two 
very different techniques. For maximum sensitivity to sinu- 
soidal pulse profiles, we performed a coherent Fourier analysis 
of all the zero-order events from the complete DDT observa- 
tion. For better sensitivity to more complicated pulse profiles, 
we performed period folding searches on the much smaller set 
of time-filtered events with the significance of each trial deter- 
mined by the Bayesian method developed by Gregory & Loredo 
(1992, 1996). 

3.1. Fourier Analysis 

We Fast Fourier Transformed the 420 million point time se- 
ries described in §2 and searched the resulting Fourier ampli- 
tudes using an advanced pulsar search code. The search in- 
cluded harmonic summing to improve sensitivity to low duty- 
cycle pulsations, Fourier interpolation to minimize the effects 
of "scalloping" (van der Klis 1989), and the ability to compen- 
sate for signals with a constant frequency derivative (i.e. an "ac- 
celeration" search) by matched filtering of the complex Fourier 
amplitudes with a series of template responses (Ransom et al. 
2001; Ransom 2001). 

Due to the limitations in the HRC-S time resolution, we lim- 
ited our search to frequencies / < 100 Hz. Similarly, we re- 
stricted the range of acceptable frequency derivatives to be from 
-5 x 10~ 10 < / < OHzs" 1 , which encompasses the range ob- 
served for all known isolated pulsars. Assuming spin down via 

5 See http://asc. harvard . edu/ cal /Links /Hrc/CIP /timing . html 



magnetic dipole radiation, pulsars with period P and magnetic 
field strengths B > 3.2 x 10 19 P 3 / 2 r _1 Gauss will drift across 
multiple Fourier bins during observations of duration T sec- 
onds. Un-accelerated searches are considerably less sensitive 
to pulsars with magnetic fields higher than this threshold. For 
a 100 ms period and a 630 ks time series (as for the DDT 
observation), this implies reduced sensitivity for pulsars with 
B > 1.6 x 10 12 G, thereby effectively eliminating a significant 
fraction of possible phase space for such sources. A near worst- 
case scenario for un-accelerated searches is the Crab pulsar (for 
which / = 29.82 Hz and/ = -3.7 x lO^Hzs" 1 ) which would 
drift by |/T 2 | ^149 Fourier bins during the observation. In 
a raw power spectrum such a signal would be smeared below 
detectability. After accounting for the number of independent 
trials searched, no candidates were detected with an equivalent 
Gaussian significance of greater than 2 a. For completeness, 
we also searched for signals with positive frequency derivative 
from < / < 5 x 10" 10 Hz s" 1 and detected no significant pul- 
sations. 

In addition to searching the Chandra DDT observation, we 
searched the archival LETGS observation (Obsld 1 13, ~ 55 ks) 
of RXJ1856 as well as the 1997 October 9 ROSAT High Res- 
olution Imager (HRI) observation (~ 29 ks) of RXJ1856 using 
virtually identical techniques. Not surprisingly, no candidates 
were found with significance greater than 2 a in either case. In 
addition, none of the lower significance candidates from these 
searches matched any of the low-significance candidates from 
the DDT observation. 

3.2. Period Folding Analysis 

We conducted brute-force period folding of the time-filtered 
events from each of the individual observations using a mod- 
ified version of the Gregory & Loredo (1992) technique that 
allows searching over / as well as /. This method maintains 
sensitivity to a wide variety of pulsed signals by matching the 
complexity of a signal's pulse shape with an appropriate num- 
ber of bins in the folded profile aligned at the optimal phase. 
Successful use of this search algorithm has been made by Za- 
vlin et al. (2000) and Hambaryan et al. (2002). 

Each individual observation was searched using the time- 
filtered events described in §2 over a range of frequencies 
from 0.001 < / < 100 Hz and frequency derivatives from -5 x 
jq-10 < y < OHzs" 1 . The candidate lists from each observa- 
tion were then compared to those from each of the other ob- 
servations to find likely matches in / and /. No interesting 
candidates were found. 

Most of the CPU time used in the searches was spent search- 
ing the two longest observations (Obslds 3380 and 3381) to- 
gether over a range of frequencies 0.001 < / < 500 Hz and fre- 
quency derivatives from -5 x 10~ 10 < / < Hz s" 1 . The num- 
ber of operations required for a folding search that includes / 
goes as ~ NphotT 2 , where N p hoi is the total number of photons 
and T is the total duration of the observation (N p hot = 3582 and 
T = 395 ks for these two observations together). This search 
required over 10 16 operations and produced no statistically sig- 
nificant candidates. 

3.3. Pulsed Fraction Limits 

Since no candidates were found in any part of our search, 
upper limits on the pulsed fraction of RXJ1856 can be derived 
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based on the predicted response of a worst-case sinusoidal sig- 
nal during the Fourier analysis. We use the standard defini- 
tion of the pulsed fraction f p = a/(a + b), where a and b are 
the pulsed and unpulsed count rates respectively. Since a si- 
nusoidal signal in the presence of noise with a total number of 
events N p h t =(a+b)T produces an expected normalized power 6 
of (P) = fpN p ho,/4-+ 1, the pulsed fraction can be written as 

f p =[4((P)-l) /N phot ] 1/2 (see e.g. van der Klis 1989; Vaughan 
et al. 1994). 

When no pulsation is detected during a search, an upper limit 
on the pulsed fraction that the data could still contain can be cal- 
culated at some level of confidence C, based on the maximum 
observed power in the search. Vaughan et al. (1994) describe 
how to calculate the signal power that would be required to pro- 
duce a measured power greater than the maximum observed 
power a fraction C of the time. These calculations are sensi- 
tive to how the data are binned and the number of independent 
trials searched — quantities that can be difficult to estimate. 
Once this signal power has been calculated, substituting it for 
((P)-l) in the definition for pulsed fraction gives the upper 
limit on f p at a confidence level C. 

Our best estimates for the limiting pulsed fraction for 
RXJ1856 based on the new Chandra observations, as well as 
the archival Chandra and ROSAT observations, are given in Ta- 
ble 2. Slightly different upper limits are quoted based on the 
fact that the number of trials is very different when compar- 
ing an acceleration search with a more standard un-accelerated 
Fourier analysis. For the archival data sets, our 50% confidence 
upper limits for un-accelerated searches are roughly consistent 
with (but slightly more conservative than) the values reported 
in Burwitz et al. (2001) and Pons et al. (2002). For "normal" 
pulsars with relatively slow spin periods (i.e. / < 5 Hz and 
/ ^ 10" 12 Hzs _1 ), the un-accelerated values are more appro- 
priate. 

4. DISCUSSION 



The 99% confidence upper limits determined here of 4.1% 
or 4.5% for un-accelerated or accelerated signals respectively, 
are the most constraining limits yet on the pulsed fraction of X- 
ray emission from RXJ1856. The fact that no pulsations have 
been observed in such sensitive observations is a surprising re- 
sult, especially considering the fact that three other isolated NS 
candidates have measured pulsed fractions of well over 10% 
(see §1). There are at least three potential reasons for the lack 
of observed pulsations: 1) Radiation from the NS is uniform 
in intensity and emitted isotropically (possibly due to a uni- 
form temperature distribution and a weak magnetic field). 2) 
The viewing angle or emitting geometries are unfavorable. 3) 
Gravitational bending effects near the NS decrease the intrinsic 
pulsed fraction. We will address each of these possibilities in 
turn. 

While uniform temperature blackbody models with T ~ 
60 eV seem to provide the best fits to the observed X-ray spec- 
tra, for the ~ 60 pc distance of Walter (2001) the implied NS 
"radiation" or "redshifted" radius R^ = R/{\ - 2GM/Rc 2 ) l l 2 
is too small (R^ <~ 2 -4km) for any known EOS (Burwitz 
et al. 2001; Pons et al. 2002). While the new distance mea- 

6 This normalization is a factor of two smaller than that used in some other X-ray timing papers. 

7 The definition of pulsed fraction as given in eqn. 16 of Psaltis et al. (2000) is equivalent to that given in §3.3. 

8 Note that the axis of ordinates in this figure is mislabeled and should be a factor of ten lower. 

9 Our definition of p, using the more traditional units km Mri is equivalent to that of Psaltis et al. (2000), who set G = c = 1 giving p = Rc 2 /2GM. 



surement of Kaplan et al. (2002) increases the radiation radius 
to /?oo ^5-9 km, since the true radius is always less than the 
radiation radius these values are still difficult to reconcile with 
most EOSs — although self-bound quark configurations (i.e. 
strange stars) may be possible (Lattimer & Prakash 2001). Un- 
fortunately, an even bigger problem with the uniform tempera- 
ture blackbody models is that they under-predict the measured 
optical flux by a factor of ~ 2 (Pons et al. 2002). 

If the surface temperature distribution of a NS is non- 
uniform, we may expect to see X-ray pulsations at the star's 
rotation period resulting from the motion of the hotter surface 
features (presumably polar cap "hotspots") through our line-of- 
sight. There are at least two geometries, though, that preclude 
the observation of such pulsations no matter what the proper- 
ties of the NS or its hotspots: when the line-of-sight is aligned 
with the rotation axis of the star, or when the polar cap axis 
is aligned with the rotation axis. Neither of these possibilities 
can be ruled out with the current observations, and given the 
modulation strengths of the other observed isolated NSs, this 
may be the simplest explanation for the lack of pulsations from 
RXJ1856. 

Finally, gravitational bending effects on pulsations from 
hotspots on neutron stars have been studied in detail by many 
groups (e.g. Pechenick, Ftaclas, & Cohen 1983; Riffert & 
Meszaros 1988; Zavlin, Shibanov, & Pavlov 1995; Psaltis, 
Ozel, & DeDeo 2000). All have showed that the measured 
pulsed fraction is a strong function of the angle between the 
rotation axis and polar cap axis, the angle between the rotation 
axis and the line-of-sight, the hotspot size as a fraction of the 
total surface area (a), and the mass and radius of the neutron 
star itself. 

For RXJ1856, one can compensate for the under-prediction 
of the optical/UV flux by the uniform temperature blackbody 
model (when fit to the X-ray data alone) with the simple addi- 
tion of a cool (~ 20 eV) blackbody component from the entire 
NS surface. In this model, the X-ray flux comes from ther- 
mal emission from much smaller (and presumably polar cap) 
hotspots (Burwitz et al. 2001; Pons et al. 2002). For the 60 pc 
distance to RXJ1856, Pons et al. (2002) find acceptable two- 
component blackbody fits to the optical, UV, and X-ray data 
with a ~ 0.2, provided that R^ < 10km and M < 1 .3 M Q (due 
to the "causality" limit). If such stars are possible, those with 
"reasonable" masses (i.e. M > 1.2M Q ) would be extremely 
compact, with R/M < 6kmMg. The extreme gravitational 
fields of such compact NSs significantly bend the light emit- 
ted from the surface and reduce the observed pulsed fraction. 

Psaltis et al. (2000) calculated the likelihood of detecting pul- 
sations of a given pulsed fraction 7 after integrating over all pos- 
sible polar cap and line-of-sight angles as a function of a and 
the neutron star radius-to-mass ratio. Figure 6 8 of Psaltis et al. 
(2000) shows the fraction of neutron stars with pulsed fraction 
greater than a specific value, plotted for several values of the 
radius-to-mass ratio 9 p <~ 3/J/MkmMQ 1 . For stars as com- 
pact as those implied by the 60 pc two-component fit mentioned 
above, less than <~ 20% yield pulsed fractions at or above the 
limits reported here. 

By transforming the two-component blackbody fit to 140 pc 



and noting that a is independent of the distance, we find R^ < 
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22 km. This upper limit is compatible with virtually all modern 
EOSs which predict 12 < R^ < 16 km (Lattimer & Prakash 
2001). For a "canonical" neutron star with R ~ 11.5 km, M = 
1.35 M Q , and R/M ~ 3p ~ gkmlVlg 1 , the majority (~ 80%) 
would show pulsed fractions in excess of our limits (Psaltis 
et al. 2000). Similarly, except for highly unlikely and unfa- 
vorable viewing geometries, stars with much higher values of 
R/M (i.e. > 12km Mq 1 ) may be excluded by the lack of pulsa- 
tions from RXJ1856. These facts imply — but do not require 
— that RXJ1856 is a relatively compact NS. 

It is important to remember in this discussion of NS hotspots 
that no truly reliable estimates of a NS radius to mass ratio can 
be derived from blackbody spectral fits, since these fits only 
yield estimates of the radiation radius (/?oo) provided the dis- 
tance is known. Furthermore, the inferred area of the X-ray 
emitting region is strongly dependent upon the presence and 
properties of any atmospheric component. Since there are as yet 
no atmospheric models consistent with the current data (Bur- 
witz et al. 2001), clearly, detailed spectral analysis and mod- 



eling of the 2001 -epoch Chandra data should be of great use 
in trying to explain the lack of pulsations from and the overall 
nature of this enigmatic object. 
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Table 1 

Chandra LETG/HRC-S OBSERVATION LOG 





?START 


Exposure 


Events 


Events 


Obsld 


(MJD) 


(ks) 


(Raw) 3 


(lmsFilter) b 


113 


51613.33 


55.5 


10436 


536 


3382 


52190.36 


102.0 


20072 


1094 


3380 


52192.22 


167.5 


33703 


1743 


3381 


52194.81 


171.1 


34571 


1839 


3399 


52197.50 


9.3 


1788 


104 



a These are the unfiltered events extracted from a ~ 2" re- 
gion around the nominal position of RXJ1856: RA(J2000) = 
18 h 56 m 35. s 5, DEC(J2000) = -37° 54' 36f'8 

b These events were time-filtered such that their arrival 
times differed by no more than 1 ms from the event recorded 
after them (see §2). 



Table 2 

Pulsed Fraction Limits for RX J1856. 5-3754 







No / Search 






With /' Search 




Obslds 


p 

1 max 


Ntrials 


Jp 


f99% 
Jp 


p 

1 max 


^ trials 


Jp 


f99% 
Jp 


400864 3 


18.4 


2.8 


< 7.0% 


< 9.8% 
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Note. — These are the pulsed fraction upper limits for RX J 1856.5-3754 based on the Fourier 
searches discussed in §3.1. The highest normalized power found in each search is denoted by P max 
and the approximate number of independent trials (including / trials and gaps in the data) is given by 
Ntrials in units of 10 6 trials. 



